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ABSTRACT 

We derive the density profile for coUisionless dissipationless dark matter haloes in hi- 
erarchical cosmologies making use of the Secondary Infall (SI) model. The novelties 
are: i) we deal with triaxial virialised objects; ii) their seeds in the linear regime are 
peaks endowed with unconvolved spherically averaged density profiles according to the 
peak formalism; iii) the initial peculiar velocities are taken into account; and iv) ac- 
creting haloes are assumed to develop from the inside out, keeping the instantaneous 
inner system unaltered. The validity of this latter assumption is accurately checked by 
comparing analytical predictions on such a growth with the results of numerical sim- 
ulation. We show that the spherically averaged density profile of virialised objects can 
be inferred with no need to specify their shape. The typical spherically averaged halo 
density profile is inferred, down to arbitrarily small radii, from the power-spectrum 
of density perturbations. The predicted profile in the ACDM cosmology is approx- 
imately described by an Einasto profile, meaning that it does not have a cusp but 
rather a core, where the inner slope slowly converges to zero. Down to one hundredth 
the total radius, the profile has the right NFW and Einasto forms, being close to the 
latter down to a radius of about four orders of magnitude less. The inner consistency 
of the model implies that the density profiles of haloes harbour no information on 
their past aggregation history. This would explain why major mergers do not alter 
the typical density profile of virialised objects formed by SI and do not invalidate the 
peak formalism based on such a formation. 
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1 INTRODUCTION 

In the last two decades, observations and, particularly sim- 
ulations, have provided us with detailed information on the 
structure and kinematics of bottom- up hierarchically assem- 
bled virialised dark matter haloes. However, from the theo- 
retical viewpoint, the situation is far from satisfactory. The 
way all these properties settle down remains to be elucidated 
and their connection with the power-spectrum of density 
perturbations is unknown. 

The accurate modelling of virialised self-gravitating col- 
lisionless dissipationless systems is an old unresolved is- 
sue. Most efforts have focused on determining the equi- 
librium density profile for objects formed from the mono- 
lithic collapse of an isolated spherically symmetric seed 
with outward-decreasing density profile and pure Hubble 
velocity field, the so-called Sec ondary Infall ( SI) m odel. 
Following the seminal work by iGunn fc GottI l|l972l 'l. the 
power-law density profile was derived under the self- 
similar approximation and/or making use of an adiabatic 
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invar iant during virialisation, both for pure radial or- 
bits (Gunn"l977: 'Fillmore fc Gol dreichI [ l984l : iBertschingeij 
1985; Hoffman^ fc Sh aham 198?^) and non-radial orbits 
Ryden fc Gunnll 19871: Iwhite fc Zaritskvlll992l : lNusseill200ll : 



iLe Delliou fc Henriksenll2003h . The departures from spher- 
ical symmetry (by adopting, instead, cylindrical symme- 
try; fRvdcn'T99d; Moutardc ct al. 1995) and self-similarity 



try; iKyacn i 

(|Lokas 2000 : .Lokas fc Hoffm an 2000) were also investigated. 

These theoretical results were tested and comple- 
mented by the information draw n from specifically de- 
signed numerical tools Gottj (|l975f) : IWiUiams ct al. (200|) 
as well as full cosmological A^-body simulations (e.g . 
Frenk et al.lll985l: lOuinn et al.|[T98^ : lEfstathiou et al.lll988l : 



Zaroubi et al.lll99d ). One important finding along this lat- 



ter line was that cold dark matter (CDM) haloes with dif- 
ferent masse s show similar scaled spher i cally averaged den - 
sity profiles JPubinski fc Carlberg|[l99ll : ICrone et al.lll994l ') . 
iNavarro et al.l ( 19971 ) showed that they are well- fitted down 
to about one hundredth of the total radius by a simple an- 
alytic expression, the so-called NFW profile, that deviates 
from a power-law. This finding opened a lively debate about 
the value of the central asymptotic behaviour of the halo 
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199S 



Moo re et al 

Ghigna et al. 2000 : Jing fc Suto 2000 ; Power et al 



densi t y profile (e.g . Fukushige 



200a : lHavashi et al.ll2004 l. More recently. the lEinastol jl965l ) 



profil e was shown to give even better fits down to small er 
radii (|Navarro et al.ll2004 [ioiol : iMerritt et"aLll2005l . [200^ ). 

The origin of this density profile remains to be un- 
derstood. Certainly, it can only arise from the way haloes 
aggregate their mass, which, in hierarchical cosmologies, 
is through continuous mergers. But the dynamical ef- 
fect of mergers depends on the relative mass of the cap- 
tured and capturing haloes, AM/M. For this reason, it 
is usually distinguished between minor and major mergers 
(AM /M greater or less than about 0.3: ISalvador-Sole et al.l 
Major mergers have a dramatic effect each on 
the structure of the object. While minor mergers con- 
tribute jointly, together with the capture of diffuse mat- 
ter (if any), to the so-called accretion that yields only a 
small secular effect on the accreting object. Some authors 
(Svcr fc Wh"it3 1199 8': Raig et al] 1 19981 : Subramanian et "al] 
[2000 ; Dekel et al. 2003) studied the possibility that the non- 



of repeated major mergers. Others ( 


Avila-Reese et al. 19981: 


Huss et alj 19991: Del Popolo et al 


2OOOI: lAscasibar et all 


2004 iMacMillan et al.|l2006'") concentrated on the effects of 



pure accretion (PA), that proceeds according to the sim- 
ple SI model above from peaks (secondary m axima) in the 
primordial random Gaus sian density field (|Doroshkevichl 
ll97d : lBardeen et aLllTgsi ). The density profile found in this 
latter scenario appears to be similar, indeed, to that of 
simulated haloes. Furthermore, recent cosmological simula- 
tions have confirmed that major mergers play no central 
role in setting the structure and kinematics of dark matter 
haloes (|Wang fc W hite 2009). However, the reason why ma- 
jor mergers would not alter the density profile set by PA is 
not understood. 

In the present paper, we develop an accurate model of 
halo density profile that clarifies all these issues. For sim- 
plicity, we consider pure dark matter systems, that is we 
neglect the effects of baryons on halo structure. This will al- 
low us to directly compare the theoretical predictions of the 
model to the results of A'^-body simulations. This model is 
built within the SI framework and uses the peak formalism, 
which assumes there is a one-to-one mapping between haloes 
and the density peaks. However, peaks are not assumed to be 
spherically symmetric and e ndowed with the i r typi cal filtered 
density profile calculated bv lBardeen et al~ 1 l|l986l . hereafter 
the BBKS profile ) as usual, but we take into account that 
they are triaxial (|Doroshkevichl Il970l : iBardeen et al.l 1 19861 ) 
and consider their accurate spherically averaged unconvolved 
density profile. In addition, we account for the initial pe- 
culiar velocities believed to affect the central halo density 
profi le (|Hiotelisll200^ : lAscasibar et al.|[200i : iMikheeva et all 
I2OO7I ). Finally, instead of making use of the usual adiabatic 
invariant, poorly justified near turnaround and greatly com- 
plicating the analytic derivation of the final density profile, 
we take advantage of the fact that accreting haloes develop 
outwardly, keeping their instantaneous inner structure un- 
altered. The validity of such an assumption, hereafter sim- 
ply referred to as inside-out growth, is accurately checked 
by verifying that the trends observed in simulated haloes 
evolving by smooth accretion are reproduced by an analytic 
model where the haloes are forced to grow inside-out. 



And what about the fact that major mergers are ig- 
nored in SI? This important caveat affects not only the 
modelling through SI of the halo density profile, as men- 
tioned, bu t also of halo statistics i n the so-called peak 
formalism jPeacock fc Heavens! [T990l: lBond fc Mverslll99ll : 



iManriaue fc Salvador-Solelll995l. Il 996'. hereafter MSSa and 
MSSb, respectively; iManrigue et al.. 1998.) which relies on 
the Ansatz that there is a one-to-one correspondence be- 
tween haloes and peaks as in spherical coUapscQ- Although 
no such one-to-one corres pondence is actually found in simu- 
lations (|Katz et al.lll99j ). this is because peaks show nested 
configurations (MSSa) which are not corrected for. If one 
concentrates in simply checking whether or not halo seeds, 
hereafter also called protohaloes, coincide with peaks, then 
the answer is positive in all but a few percent of cases 
compatible with the frequency of non-fuUy virialised haloes 
IPorciani et a]||2002l 'l. However, there is still the problem 
that, in principle, major mergers should blur the correspon- 
dence between peaks and haloes arising from SI. In the 
present paper, we argue that the reason why halo statis- 
tics and the spherically average density profile of virialised 
haloes seem to be unaffected by major mergers is that virial- 
isation is a real rela:xation process. We stress that the model 
deals with fully virialised haloes. This is important because, 
after a major merger, the density profile of a halo spends 
some time to adopt a stable density profile. 

The outline of the paper is as follows. In Section [21 we 
present some general relations holding for all systems regard- 
less of their symmetry. Taking into account these relations, 
we develop in Section [3] the model for triaxial haloes grown 
by PA. In Section [l] we calculate the properties of typical 
protohaloes and use them in Section [5] to derive the typical 
spherically averaged density profile for CDM haloes. In Sec- 
tion [6l we discuss the effects of major mergers. The main 
results are summarised in Section [T] 

Throughout the p aper, we adopt the ACDM Wmap7 
(jKomatsu et al.l I2OIII ) concordance model. In spite of 
this, we use Newtonian dynamics with null cosmologi- 
cal constant as its effects are irrelevant at the scale of 
virialised haloes. A package with the numerical codes 
used in the present paper is publicly available from 
www . am.ub . es/~cosmo/haloes&peaks . tgz. 



2 GENERAL RELATIONS FOR SPHERICAL 
AVERAGED PROFILES 

In the present section, we derive some general relations for 
spherically averaged quantities that hold for any arbitrary 
system regardless of its symmetry and that will later be used 
to build the model. 

Consider a self-gravitating system with arbitrary mass 
distribution, aggregation history and dynamical state. The 
local density and gravitational potential can be split in the 
respective spherical averages around any given point and the 
corresponding residuals. 



p(r,e,ip) = {p)(r) + Sp(r,e,v) 



(1) 



^ The excursion set formalism is also based on a one-to-one cor- 
respondence between haloes and overdense regions in the initial 
density field as in spherical collapse. 
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<E.(r,0,^) = (<E-)(r) + 5$(r,e,<^). (2) 
The spherically averaged gravitational potential, 

($)(,.) = A^j desine<J>(r,e,<^), (3) 

satisfies, by the Gauss theorem, the usual Poisson integral 
relation for spherically symmetric systems, 

d{$)(r) _ GM(r) 



dr ' 
where M{r) is the mass within r, 

M{R) = 4n drr'^{p){r). 



(4) 



(5) 



Taking into account the null spherical averages of Sp 
and 5$ (see eqs. [l]-[2]), the potential energy within the 
sphere of radius R, 

1 /'27r 

W{R)^- drr^/ d<fi de sine p{r,e,<p)${r, e,ip) , (6) 
^Jo Jo Jo 

can be rewritten as 

W{R) = 2tx I dr [(p) (r) (-J-) (r) + 5$)] 
Jo 







47r/ drr'(p)(r) ^^^^^^ +27r / drr^{Sp5^) 







>V(i?) + <5>V(i?) , 



(7) 



where we have introduced the "spherical" potential energy 
within R, 



W{R) = -47r / drr^(p)(r) 



GM{r) 



(8) 



The second equality in equation ([7]) follows from partial in- 
tegration of the first term on the right of the first equality, 
then application of the relation (Q, and one new partial 
integration, choosing the origin of the spherically averaged 
potential so as to have 

GM{R) 



{$)iR) = 



R 



(9) 



Note that, for spherically symmetric systems, this bound- 
ary condition coincides with considering the usual potential 
origin at infinity and the system truncated at R. 
The kinetic energy within R is 

K{R) = 2nf drr'^{p){r)a^{r) , (10) 

^0 

being (7^(r) the velocity variance in the infinitesimal shell at 
r0 Thus, the total energy of the sphere, E{R) = K{R) + 
W{R), takes the form 



E{R) = ^-k\ drr{p){r) 



a^{r) GM{r) 



(11) 



+ 27r/ drr'^{5p5<i>){r) = £{R) + 5£{R), 
Jo 

where we have introduced the "spherical" total energy 



^ Such a velocity variance coincides with the spherical average of 
the local value. Therefore, equation l|10|l could also be written in 
terms of the local velocity variance in angular brackets. 



£{R) = Ati I drr{p){r) 



Halo Density Profile 

s^{r) GM{r)' 



(12) 



Note that £[R) is written in terms of the "spherical" veloc- 
ity variance s'^(r), different, in general, from the ordinary 
velocity variance, 

(j^{r) = s'(r) + 5s'{r), (13) 

through the residual Ss'^{r) to be specified, contributing to 
the residual 5£{r). In other words, there is some freedom in 
the definition of s^ir). We will come back to this point later. 

So far we have considered a system with arbitrary mass 
distribution (symmetry) , aggregation history and dynamical 
state. If the system is in addition in equilibrium, then mul- 
tiplying the steady coUisionless Boltz mann equation (e.g. 
eq. [4p-2] in lBinnev fc Tremainelll987f ) by the radial parti- 
cle velocity and integrating over velocity and solid angle, we 
are led to 



d({p)a?) , {p){r) 



dr 



[3a,^(r)-a^(r)] 



1 

47r 



dip! de sine p{r,e,(p)dr^{r,e,(p) , (14) 



where dr stands for radial partial derivative[f|. Taking into 
account equation (|4]), equation (|14ll adopts the form 



dr r 

= ^{p){r)^^-{SpdriS^))(r), 



(15) 



identical to the Jeans equation for spherically symmetric 
systems in equilibrium but for the last term on the right. 
Multiplying equation (llSp by 47rr'^ and integrating over the 
sphere of radius R, the same steps leading to the scalar virial 
relation for spherically symmetric self-gravitating systems 
now lead to 



4TTR^{p){R)a^{R) - 2K 



-4:71 1 drr (p) (r) 
'o 



GM{r) 



-4n drr^{5pdrS<i>){r). (16) 



(17) 



Therefore, defining the so-called "spherical" radial velocity 
variance s^{r) through 

o-j (r) = Sr(r) + Ss^{r) , 



with 

5sr(r-) = 



dfr^[Ss'^{r)-f{5pdr5$){f)], (18) 



the virial relation (|16p takes the usual form for spherically 
symmetric systems, 

4-kR^{p){R)s'^,{R) 



4n I drr {p){r) 
lo 



2 GM{r) 
s (r) 



(19) 



from now on called the "spherical" virial relation. Further- 
more, defining the "spherical" scaled surface pressure term 
(from now on simply spherical surface term), S{R), equal to 

^ To derive equation I I14II it is only needed such conventional 
conditions as the continuity of the local density and mean velocity 
and the fact that the velocity distribution function vanishes for 
large velocities. 
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the member on the left of equation (|19|) over the absolute 
value of W{TZ), the spherical virial relation (|19[l adopts the 
usual compact form 



2£{R) 
W{R) 



S{R). 



(20) 



As can be seen, all the ordinary quantities X can be ex- 
pressed in the form X = X + SX with the quantities X hav- 
ing the same form for spherically symmetric systems, hence 
why they are labelled "spherical" . The residuals SX always 
measure the deviation from sphericity of the quantities X 
in the sense that, when 5p and 5<E> vanish, 5X also vanish 
and X become equal to the spherical counterparts X and, 
hence, recover the form for spherically symmetric systems. 
The dimensionless quantities \SX/X\ are, however, not nec- 
essarily less than one, so we should not look at the spherical 
quantities X as Oth-order approximations of X in expansion 
series on small deviations from sphericity; they are fully ex- 
act quantities and so are also the relations between them. In 
particular, equations ([5]), (I12p and (|20|l respectively giving 
the mass, spherical total energy and spherical virial relation 
or, equivalently, the relations 



{p}{r) = 



1 dM 



s^(r) = 2 



and 



d£/dr _^ GM{r) 



dM/dr 



(21) 



(22) 



(23) 



2 _ 2g(r)-W(r) 
' r dM/dr ' 

following from their differentiation are exact and hold un- 
changed regardless of the particular shape of the object. 
Thus, the profiles {p){r), s'^{r) and Sr(r) do not depend on 
shape and can be determined by solving equations (PTjl - lP^ . 
In contrast, the ordinary kinematic profiles (T^{r) and Ciir) 
cannot be obtained from a similar set of equations because 
they depend explicitly on the shape of the object. 

We want to remark that, even if the spherical quantities 
at any given moment behave as their ordinary counterparts 
in spherically symmetric systems, the relation between their 
values in different epochs may not be the same as in such sys- 
tems. For instance, the spherical total energy inside spheres 
of fixed mass may not be conserved in the absence of shell- 
crossing]. Only one specific choice of s^(r) will ensure such 
a conservation. Note also that, for any arbitrary mass dis- 
tribution, it might not be possible to have one only profile 
s^(r) satisfying that condition for the whole series of embed- 
ded spheres the system decomposes into. But, in the case of 
a centred triaxial system, this is always possible thanks to 
the ordered radial mapping of the mass distribution. From 
now on, we assume such a symmetry and the conservation 
of the spherical total energy in the absence of shell-crossing. 
This is crucial for the model as it will allow us to fix the 
density profile for a virialised object from the properties of 
its seed (see Sec. [Sjl . 



* If one chooses s^(r) = (j^(r), the spherical total energy is not 
conserved in ellipsoidal systems because any given sphere ex- 
changes energy with the rest of the system. 



3 THE MODEL 

We will concentrate, in a first step, on triaxial systems evolv- 
ing by PA. This includes the case of accretion along fila- 
ments. 



3.1 Inside-out Growth 

The model relies on the assumption that accreting objects 
grow inside-ou t. This is consistent with the properties of sim- 
ulate d haloes ([Salvador-Sole et al ]|l998l:ISalvador -Sole et al.l 



20051: iHelmi et all l2003l: iRomano-Dfaz et all l20o'i 



Wang et al 



2011 



In iManrigue et al.l 
(|200'iD . it was shown 



1 20031 ) and 



Salvador-Sole et all (|2007D . it was shown that the as- 
sumption that (spherical) haloes grow inside-out at the 
typical cosmological accretion rate automatically leads to a 
density profile a la NFW. Nonetheless, to be fully confident 
about it, we show below that such a growth is naturally 
expected and in agreement with the detailed results of 
simulations. 

In SI, the isodensity contours in the triaxial seed 
expand, without crossing each other, until they reach 
turnaround. This expansion is achieved in linear regime with 
all particles moving radially, so the axial ratios of shells 
are conserved. However, after reaching turnaround, parti- 
cles collapse and rebound non-radially, so the shape of the 
isodensity contours changes. Shell-crossing causes a secular 
energy transfer between shells so that particles initially ly- 
ing in a turnaround ellipsoid reach their new apocentre at 
a smaller radius and not simultaneously. However, the en- 
ergy lost by particles in one orbit is small, so the typical 
time of apocentre variation is substantially greater than the 
typical orbital period, implying that particles belonging to 
one turnaround ellipsoid will still trace an effective apoc- 
entre surface, which for symmetry reasons will also be an 
ellipsoid although smaller in size and with different axial ra- 
tios. Repeating the same reasoning from each new apocen- 
tre ellipsoid, we have that any turnaround ellipsoid evolves 
through a continuous series of ellipsoidal apocentre surfaces 
that progressively shrink and change their shape until they 
stabilise. 

During such an evolution, apocentre ellipsoids never 
cross each other. If two of them had a common point, par- 
ticles belonging to the two surfaces at that point (with null 
velocities) would follow the same evolving orbits, so the two 
surfaces would always be in contact, but this is meaningless 
because different turnaround ellipsoids come from different 
isodensity contours in the seed which do not intersect. Thus, 
during virialisation, the system contracts orderly, without 
apocentre- crossing and when a new apocentre ellipsoid stops 
shrinking, it necessarily places itself beyond all previously 
stabilised apocentre ellipsoids. Thus, the central steady ob- 
ject necessarily grows from the inside out according to the 
gradual deposition of particles with increasingly larger apoc- 
entre surfaces. 

In the above reasoning, we assumed that orbits even- 
tually stabilise. However, new shells are constantly arriv- 
ing and crossing the growing halo, consequently, orbits only 
approximately stabilise. It is therefore worthwhile verifying 
the validity of this approximation. This can be accurately 
checked by means of the detailed follo w-up of the d ensity 
profile for CDM haloes in simulations. IZhao et al.l (|2009| ) 
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Figure 1. Tracks followed in the Mg-r^ plane by typical accreting 
haloes developing inside-out, with current masses, M, equal to 
10", 10^2, 10", iQi* and lO^^^ Mq (from bottom to top) at 2 = 0, 
2, 4, 6 and 8 (circles essentially from right to left, respectively in 
black, green, red, blue and magenta). The straight dashed red line 
shows the uni versal direct i on fol lowed by such tracks, identical to 
that found bv lZhao et al in full cosmological simulations. 

The black straight line shows t he Mb—Tb rel atio n obta ined at 2 = 
also identical to that found in lZhao et al.l al. l|2009l : see Fig. 22). 



and iMunoz-Cuartas et al] (|201lh carried out such a follow- 
up and found that the rs and Ms NFW shape parameter^ 
of accreting haloes are not constant but vary slightly with 
time. This was interpreted as evidence that the inner struc- 
ture of accreting haloes is changing. But this is not the only 
possible interpretation. As the densityprofile of haloes is 
not perfectly fitted by the NFW profilqj, even if the inner 
density profile did not change, its fit by the NFW law out to 
progressively larger radii should result in slightly different 
best values of Vs and Ms. 

To confirm that our interpretation for this change is 
correct we have followed in the Mb-Tb log-log plane the evo- 
lution of accreting halos forced to grow inside-out at the 
typical accretion rate give n by the excur sion set formal- 
ism according to. Salvador- Sole et al.l ()2007l) model. The re - 
sult, in the same cosmology as used in IZhao et all l|2009l l. 
is shown in Figure [1] Even if, by construction, haloes grow 
inside-out without changing their instantaneous inner struc- 
ture, when their density profile is fitted by a NFW profile, 
they are found to move in the Mg-rg log-log plane along 
one fixed direction over a distance that increases with cur- 
rent halo mass. This behaviour is identical to that found 



° Mb is defined as the mass within the scale radius rg where the 
logarithmic slope of the density profile is equal to —2. 
^ The larger the mass, th e closer the density pro file is to a power- 
law with index ~ —1.5-2 jTasitsiomi et al.|[2004) . 



bv lZhao et al.l : accreting halos moved in the Mg-rg log-log 
plane along straight lines with exactly the same universal 
slope {Mb oc rg'^^) as in our experiment and ended up at 
z = lying along another straight line, with a different 



slope (Ms 



) also identical to that found in our ex- 



periment (see their Fig. 22). The fact that our constrained 
model reproduces the observed trend convincingly demon- 
strates haloes primarily grow inside-out. 

3.2 Radius Encompassing a Given Mass 

As we will show, by assuming that haloes grow inside-out 
via PA and that the spherical total energy in spheres of fixed 
mass is conserved in the absence of shell-crossing, we are able 
to infer the radius enclosing a given mass in a halo from the 
spherical energy distribution of its progenitor protohalo. 

To do this, we will consider the virial relation (|20[l 
for a perfectly uniform sphere with mass M, which implies 
W(M) = -3GM^/{5R), with spherical total energy £{M) 
equal to that of the system at turnaround £ta.{M) and null 
spherical surface term S{M), 



R{M) 



3 GM^ 



10 SUM) 



(24) 



This equation is the often used estimate for the radius en- 
compassing mass M in spherical systems with an unknown 
internal mass distribution (Bryan and Norman 1998). In 
principle, R{M) is not believed to give an exact measure 
of that mass in the real virialised object as this has non- 
uniform density profile, its spherical total energy is equal to 
S{M) instead of £ta.{M) and its surface term is not null. Yet, 
as shown below, as a consequence of the inside-out growth 
of accreting objects formed by PA, the inaccuracies above 
exactly cancel and both radii turn out to fully coincide. 

To see this, we will deform the system since its shells 
reach turnaround so as to construct a virialised toy ob- 
ject satisfying the conditions leading to equation (|24p . 
This is only possible in PA where the equivalent radius 
of turnaround ellipsoids increases with increasing time and 
the central virialised object grows inside-out. Shells reach- 
ing turnaround can then be virtually contracted one after 
the other without any crossing so as to match the mass pro- 
file M{r) (although not necessarily the ellipsoidal isodensity 
contours) of the real virialised object developed until some 
time t. By "virtual" motion we mean a motion of shells that 
preserves their particle energy and angular momentum, but 
is disconnected from the real timing of the system. Of course, 
such a motion will not recover at the same time the axial ra- 
tios of the isodensity contours of the virialised object, but we 
only need to recover the mass profile, which by conveniently 
contracting each new shell is guaranteed (there is one degree 
of freedom: the final equivalent radius of the contracted el- 
lipsoidal shell and one quantity to match: the mass within 
the new infinitesimally larger radius). By construction, the 
new toy object so built has the same mass profile M{r) (al- 
though not the total energy profile E{r) due to the different 
ellipsoidal mass distribution) as the virialised object. But 
the spherical total energy in such a toy object is equal to 
that at turn-around, £ta.[M{r)], as there has been no shell- 
crossing and £ is conserved for the appropriate choice of 
s^(r) (see SecfH). 

Of course, such a toy object is not in equilibrium. But 
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the quantity fta(-R) — with W{R) equal to the po- 

tential energy of a homogeneous sphere with mass Af , is 
positive (see eq. [25] below). Therefore, we can virtually ex- 
pand each ellipsoidal isodensity contour of the toy object, 
from the centre out to the edge avoiding shell-crossing, so 
as to end up with a uniform density equal to the mean den- 
sity p{R) of the real object inside R and still have an excess 
of spherical kinetic energy. Thus, this can be re-distributed 
over the sphere and the radial and tangential components 
of the spherical velocity variance, s{r), can be locally ex- 
changed so as to satisfy the spherical virial relation (|19|) at 
every radius. This is again possible to achieve because there 
are two degrees of freedom, s(r) and Sr(r), and two condi- 
tions to fulfil, the spherical total kinetic energy in excess 
and the spherical virial equation. In this way, we have built 
a steady homogeneous toy object with radius, mass and to- 
tal energy respectively equal to R, M and £ta(M) and with 
null value of Si{R) (there is no density outside the sphere 
because shells having not reached turnaround at t have not 
been contracted). Hence, this uniform toy object satisfies 
the spherical virial relation 



2SUR) 
W{R) 



(25) 



which, given the equality W(7?) = -3GA/^/(5i?), implies 
equation (|24|l . 

The exact equation (|24p allows one to determine the 
spherically averaged density profile of the virialised object 
from the properties of its seed. To do this we must simply 
take into account that, during the initial expanding phase 
of shells, there is no shell-crossing, so the spherical total en- 
ergy in spheres with fixed mass is conserved. We can then 
replace, in equation (|24|l . £ts,{M) by its value in the pro- 
tohalo, £p{M), and obtain, by inversion of r{M), the mass 
profile M{r) and, through equation (|2ip . the spherically av- 
eraged density profile (p) (r) of the virialised object from the 
spherical total energy distribution in the seed. 

Furthermore, as the functions {p)(r), s^(r) and s^{r) do 
not depend on the shape of the object, we can think in the 
spherically symmetric case to infer the latter two functions 
from the former one. In such objects, orbits are purely ra- 
dial because they collapse and virialise radially. We therefore 
have s'^(r) = s'^{r). Equations (|22|) and (|23p then lead to a 
differential equation for £{r) that can be readily integrated 
for the boundary condition £ — at r = 0, the result being 



£{R) 



-R I dr 



47r (p)(r-) GM{r) + 



W{r) 
2r2 



(26) 



Once {p){r) is known, we can calculate W(r) (eq. [6] and 
then £{R) (eq. [26]), and apply any of the two relations (|22p 
and (|23p to infer the profile s^(r). But the interest of this 
expression is only formal because, as mentioned, the profile 
s'^(r) is not observable, not even in the spherically symmet- 
ric case where s(r) = a{r) and Sr(r) = fr(r). The reason 
for this is that objects formed by PA are never spherically 
symmetric. The tidal field of surrounding density fluctua- 
tions causes sp herically symme tric seeds to undergo ellip- 
soidal collapse ijZeldovich Il970l ). Even if such a tidal fleld 
is artificially removed, the system suffers radial orbit insta- 
bility , also leading to triaxia l virialised objects (|Huss et al.l 
ll999l : lMacMillan et al.|[200^ ) . 



The expression (|26[) is nonetheless useful to compute 
the spherical total energy dissipation factor since the time 
of the protohalo, D(Af) = £{M)/£p{M), with £p{M) given 
by equation (|24|) . This has the following interesting implica- 
tion. In objects formed by PA, there is always some spher- 
ical energy loss through shell-crossing during virialisation, 
so T>{M) is greater than one. Taking into account equations 
(1201) and (|25|) . this dissipation factor can be written in the 
form 



V{M) = Y^SMl\i^s{M)] . 
W{M) 

Thus, if T>[M) is greater than one, the ratio 



W(Af) 



W(A/) 



1 + 



R(M) 



dr 



(27) 



(28) 



R{M) R^p'^iR) 

must also be greater than on^. By differentiation and after 
a little algebra it can be seen that ratio W{M)/W(M) is 
greater than one provided only p{r) is outward-decreasing. 
Consequently, we are led to the conclusion that the spheri- 
cally averaged mean inner density profile of any collision- 
less dissipationless virialised object formed hierarchically 
(through PA but also through major mergers; see Sec. [6} 
IS necessarily outward-decreasing. 



4 PROTOHALOES 

According to the previous results, to calculate the spheri- 
cally averaged density profile for a CDM halo grown by PA 
we only need the spherical energy distribution, £ta(A/), of 
its seed, namely a peak in the initial density field filtered at 
the mass scale of the halo. The cosmic time ti where the ini- 
tial density field has to be considered is arbitrary although 
small enough for the protohalo to be in linear regime. 

Such a spherical energy distribution is given, in the 
parametric form, by the spherical total energy of the proto- 
halo. 



£:p(i?p) = 4n / drpTp (pp)(rp) 



[Hirp - vp{rp)f , 4iRp) GM{rp) 



(29) 



in centred spheres of radii Rp encompassing the masses M, 



M{Rp)^47v drprl{pp){rp). 
Jo 



(30) 



In equations (|29p and (I30p . Hi is the Hubble constant at the 
cosmic time ti of the protohalo, {pp){rp) is its spherically 
averaged density profile, Vp{rp) is the peculiar velocity, to 
Oth-order in the deviations from sphericity, due to the grav- 
itational pull by the mass excess within the radius rp and 
(jp(i?p) is the uniform and isotropic backgroundpeculiar ve- 
locity dispersion inside the sphere of ra dius R^ 



The peculiar velocity Vp{rp) is (e.g. |Peeblej|l980h . 

^ As W{M)/VV{M) is positive (see eq. [28]), the factor 1~S(M) 
must also be positive and less than one. 

* The starting value of the spherical total energy can be com- 
puted taking Sp = cTp. 
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2G [M{rp) ~ 47Trlpi/3] 



(31) 



being pi the mean cosmic density at Bringing tliis expres- 
sion of Vp{rp) into equation (|29p and neglecting the second 
order term in the perturbation with Vp{rp) (the term with 
a^{Rp) may be greater than this owing to the contribution 
of o-n^j(ti)), we are led to 



£AR.)-l£^iRp) + M{Rp)^^ 



(32) 



where £^{Rp) stands for the spherical total energy of the 
protohalo in the case of pure Hubble flow (i.e. eq. [29] with 
null Vp and cTp). In principle, ap{Rp) has two contributions: 
one arising from the dark matter particle velocities at decou- 
pling, adiabatically evolved until the time ti, croMi^-Oy ^^'^ 
another one caused by rando m density fluctuations . Follow- 
ing the exact prescription bv lMikheeva et aL I (|2007l ) for this 
latter contribution, we arrive at 

al{Rp)^alM{ti) + ^f J ^P{k) 

= a^M(ii) +S? [a\{0) - 2a\{Rp) + a\{V2Rp)] , (33) 

where Ei is defined as y/2TvHaf, being ai the initial cosmic 
scale factoiF°l. and a^i(R{) is the spectral moment of order 
— 1 for the power-spectrum, P{k), with the j-order moment 
given by 



,(i?f) 



dfc 
27r2 



P{k)k 



2j+2 -k-' 

■' e 



(34) 



But the velocity variance due to random density fiuc- 
tuations, ap{Rp) — crnM(*i) (see eq. [33) ) turns out to be 
several orders of magnitude less than GM{Rp) / Rp (a factor 
~ 10~® in the concordance model adopted here). Conse- 
quently, it can be ne glected i n gp(M) ("eg. [29l'). Cont r ary to 
the claims made bylHiotelisI l[2002l ) . [Ascasibar et al.l ([2004[ ') 
and [Mikheeva et al.l ( 2007[ ) that such a velocity dispersion 
may affect the central density profile for CDM haloes. The 
former two authors reached this conclusion by analysing the 
effects of velocity dispersions that, unlike those used here, 
were not consistently infe rred from the random density fiuc- 
tuations. [Mikheeva et al.l used the same derivations for the 
dispersion as used here, however, they considered the fiuc- 
tuations within the whole protohalo, whereas we account 
for random density fluctuations within each sphere of ra- 
dius Rp. In including fluctuations from the whole proto- 
halo, Mikheeva et al. included velocities induced by exter- 
nal density fluctuations with masses possibly larger than 
the own mass of the sphere. This overestimates the veloc- 
ity dispersion in the sphere because external density fluc- 
tuations cause it a bulk velocity, not a velocity dispersion. 
And the bulk velocity of a sphere does not hamper its col- 
lap se, whereas the i nternal velocity dispersion does, as found 
by [Mikheeva et alj for small enough Rp . The only velocity 
dispersion that does not diminish with decreasing Rp and, 
hence, can really set a minimum halo mass is that intrinsic 



^ In equation II31I I. we have taken into account that the cosmic 
virial factor /(H) ~ fi"'^ is at very approximately equal to one. 

At small f;, this is very approximately equal to the cosmic 
growth factor. 



of dark matter particles, (J^t^{U). But this is negligible in 
CDM cosmologies. 

And what about the spherically averaged peak den- 
sity proflle, { p p)(rp) l The BBKS proflle calculated by 
iBardeen et al.[ l[l986l ) gives the typical (average) spheri- 
cally averaged density contrast profile for peaks filtered with 
a Gaussian window. Unfortunately, the convolution by a 
Gaussian window results in some information loss so the 
BBKS proflle cannot be used to infer the desired proflle. In 
addition, the average profile of purely accreting haloes with 
M at t may not coincide with the proflle arising from the 
average proflle of the corresponding seeds. However, we can 
still find the typical (in the sense below) spherically averaged 
density profile of protohaloes. 

As in PA haloes grow inside-out, the density proflle in- 
side every inner ellipsoid exactly matches that of one halo 
ancestor, also evolved by PA from a peak with its own den- 
sity contrast and scale. Hence, the typical (average) spher- 
ically averaged density proflle for purely accreting haloes 
with M at t must arise from a protohalo whose spherically 
averaged density contrast proflle, (5p)(f"p), convolved with 
a Gaussian window of radius Ri corresponding to the mass 
of any halo ancestor yields, at rp = 0, the density contrast 
5pii{Rf) of the peak evolving into that ancestor. 



-,1/2 



5pk(-Rf) = — 



TV^/^Rf- 



drpr-p((5p)(rp)exp 



2 KRtJ 



(35) 



According to the peak formalism (see Sec. [6|, the typical 
(most probable) trajectory 5pk(i?f) of peaks evolving into 
the series of halo ancestors ending in a halo with M at t is 
the solution of the differential equation (SSb) 

^ = -a;c(5pk, Ri) 02{Ri) Ri , (36) 

with the boundary condition defined by the halo, that is 
satisfying the relations 



5pk(t) = 5,{t) 



G{U) 
G{t) 



and 



Ri{M) = - 



3M 
An Pi 



1/3 



(37) 



(38) 



where G{t) is the cosmic growth factor, q = 2.75 is the ra- 
dius, in units of Ri, of the collapsing cloud with volume equal 
to Ad/pi, 5c[t{z)] = 1.93 + (5.92 - 0.472^ + O.O5460^)/(1 + 
0.0005682"^) is the critical linearly extrapolated density con- 
trast for halo formation at the redshift z (see Sec. [3| for 
those values of q and (5c[t(2)]). In equation (I36p . a;o(5pk,-Rf) 
is the inverse of the average inverse curvature x (equal to mi- 
nus the Laplacian over a^) for the distribution of curvatures 
(BBKS), 



-){Ri,Spi,)-- 



(2^) 



-1/2 



(1-^2)1/: 



da; — f{x) e 



(39) 



at peaks with (5pk and Ri, being 



3x 



+ 



erf 



1/2 



+ erf 



[2) 2 



31a;^ 



8 

+ 5 



; x" 8 
+ IT-5 



(40) 
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Figure 2. Typical <5pit— ^Jf (in physical units) peak trajectory 
of peaks at z = 100 (solid lines) giving rise to the series of halo 
ancestors evolving by PA into haloes with current masses M equal 
to 10, 1 and 0.1 times the critical mass for collapse, M* = 3.6 X 
10^^ M0 (curves from top to bottom, respectively in orange, red 
and brown) . 

and 172 (i?f) the second order spectral moment, where 7 and 
a;* are respectively defined, in terms of the spectral mo- 
ments, as a-f/{aoa2) and 7(5pk/o"o. This distribution func- 
tion is a very peaked, quite symmetric, bell-shaped function 
so that the function J?f(5pk) inverse of the 5pk(J?f) solution 
of equation (|36|) is traced by peaks with the average (es- 
sentially equal to the most probable) inverse curvature at 
each point {Spk,Rf). Note that the slope di?f/d(5pk trans- 
lates into the typical (average or most probable) accretion 
rate, dM/dt, of haloes with M{Ri) at t(5pk) evolving from 
those peaks. In Figure (2] we show the typical peak trajecto- 
ries at 2 = 100 leading to typical haloes with current masses 
equal to 10" ^M*, M* and 10M«, where M* is the critical 
mass for collapse in the concordance model (3.6 x 10^^ Mq). 

Given a typical peak trajectory, 5pk(/?f), equation (|35|) 
is a Fredholm integral equation of first kind for {<5p)(r-p). 
Through the changes y ~ and x = l/(27?f^), it takes the 
form of a two-sided Laplace integral transform, 

d^//(^/)e-^^ (41) 

- 00 

with f{y) and g{x) respectively equal to y^^'^ {&p){y^^^) and 
2y2^(2a;)-^''25pk[(2a;)-i''2], which can be solved in the stan- 
dard way. Extending x to the complex space and taking 
X = i27r^, equation adopts the form of a Fourier Trans- 
form, 

/CtO 
dy/(y)e-^^-«^ (42) 
- 00 

being .g(0 = 2^(i47rC)~='''^[(5fl(0 + where Sr{Ci 



Figure 3. Spherically averaged unconvolved density profiles for 
the seeds of typical purely accreting haloes with the same current 
masses M as in Figure|2]( same colours, now ordered from bottom 
to top), obtained by inversion of the peak trajectories shown in 
that Figure. Rp is the radius of the protoobject with mass M . The 
filtering radii Rf/Rp used in the convolutions shown in Figure |4] 
are marked with arrows. 

and (5/(C) stand respectively for the real and imaginary parts 
of (5pk(0 = <^pk[(4'''"^)^^''^]- Thus, equation p5|) can be in- 
verted by simply taking the inverse Fourier transform of 
equation (|42|) . 

As f(y) is a real function, the real and imaginary parts 
of g(£_) must be even and odd, respectively. Both conditions 
must be satisfied at the same time, so there are only two pos- 
sibilities: either (Sb(^) = '5/(C) or ^rH) ~ ~^i{0- The latter 
possibility leads to a nonphysical solution, so g{^) must be 
a pure imaginary odd function, i.e. (5_r(^) — 5/(C)- To deter- 
mine (5pk(0 we must solve the differential equation (I36|) in 
the complex space for ^, that is for its real and imaginary 
parts separately. But, as the two solutions Sr and Sj axe 
identical, we can directly solve it in the real space for 7?f 
and then take Sn = 3i = Jpk/%/2, after the conversion from 
Rf to ^. Once the solution (5pk(C) is known, we can readily 
calculate the function g{^) and take its inverse Fourier trans- 
form, which leads to the wanted protohalo density contrast 
profile, {5p){rp), after the change — y. 

The spherically averaged unconvolved protohalo density 
profiles, {pp}{rp), obtained from the typical peak trajectories 
depicted in Figure [2] are shown in Figure [3] In Figure 14] we 
plot the convolved density contrast profiles for the seeds 
of three arbitrary ancestors (with masses M{Rf) given by 
eq. for the three arbitrary filtering radii) of the halo 
with final mass M*, corresponding to the filtering radii in 
units of the protohalo radius, Ri/ Rp, marked with arrows in 
Figure O For comparison, we also plot the BBKS profiles for 
peaks with identical central density contrasts and filtering 
radii. As can be seen, each couple of curves is very similar; 
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Figure 4. Spherically averaged density contrast profiles for the 
seed at 2; = 100 of a halo with current mass equal to Af* (in red 
in Figs. [2] and [S]! convolved by a Gaussian window (solid lines) 
with the filtering radii quoted (in physical units) and marked 
with arrows in Figure |3] (same colours) , compared to the BBKS 
profiles with identical central density contrasts and filtering scales 
(dashed coloured lines). i?p"'^ are the radii of the seeds of the 
corresponding halo ancestors. 

the small deviation observed (of less than 10 % within the 
radius iip"'^ of each ancestor seed) is likely due to round-off 
errorj"!. Thus, there is no evidence that the peak leading 
to an object with the average spherically averaged density 
profile for haloes with M at t is significantly different from 
the average peak leading to such haloes. 

Interestingly, the peak trajectories converge to a finite 
value (null asymptotic slope) as Rf approaches to zero, as 
can be seen after some algebra from the null asymptotic 
limit of the right hand-side member of equation (|35|) . This 
in turn implies, from equation (|35p and the relation between 
{5pp){rp) and (Pp)(rp), that the unconvolved density profiles 
for protohaloes with different masses also converge to a finite 
value. This is consistent with the behaviour at small radii 
of the unconvolved density profiles of protohaloes shown in 
FigureO This finite limit is approached slightly more rapidly 
than in the case of typical peak trajectories. 



5 HALO DENSITY PROFILE 

Thus, to derive the typical spherically averaged density pro- 
file of a halo with M at t we must follow the following steps: 

The convolution is achieved by taking the product of the 3- 
D Fourier transform of the unconvolved density contrast profile 
and the Gaussian window. As the Fourier transform of the un- 
convolved density contrast profile has no compact support, some 
aliasing is present. 



1) solve equation (|36|) for the typical peak trajectory (5pk(-Rf ) 
leading by PA to the typical halo with M at t; 2) from such 
a peak trajectory, invert equation to find (5p)(rp) and 
from it the spherically averaged density profile of the pro- 
tohalo, {pp)(rp); 3) given the protohalo density profile, de- 
termine the spherical energy distribution £p{A4) in the pro- 
tohalo by means of equations (|30p and (|29|) : and 4) making 
use of this latter function, invert equation (|24|l to find the 
typical mass profile Af (r) of the halo and, from it (eq. |21)). 
the typical spherically averaged density profile {p){r). 

The density profiles so obtained (hereafter the theoret- 
ical profiles) for haloes with the same masses as in Fig- 
ures [2H3] are comp ared in Figure \S\ t o the corresponding 
NFW profiles with IZhao et al.1 l|2009l ) mass-concentration 
and a total halo r adius equal to the virial radius Rvii of 
iBrvan fc Normajil (|l998l ) (essentially equal to rgo). As can 
be seen, there is good agreement between each couple of 
curves down to about one hundredth -Rvir. The residuals 
Alog({p)) (prediction minus NFW profile) start being pos- 
itive at O.Of-Rvir, tend to diminish at intermediate radii 
(reaching slightly negative values for M«), then increase 
again at moderately large radii (except for 10 Af«) and fi- 
nally become negative near the halo edge. All these trends 
coincide with those shown by the same residuals (simulated 
halo mi nus NFW profil e ) for i ndividual haloes in the simula- 
tions bv lNavarro et all (|2004l : Fig. 1, left panels) except for 
the fact that the residuals of simulated haloes are notably 
larger, as expected (the theoretical profiles are supposed to 
correspond to average profiles). The only significant differ- 
ence is near the halo edge, where the residuals of the the- 
oretical profiles for low halo masses keep on decreasing for 
a longer radial range (see also Fig. (Tjl. This could be due 
to the different halo radius used in both works: the radius 
adopted in the present version of the model, arising from the 
fit to the halo mass mass function predicted in the excursion 
set formalism (see SecO is likely cl oser to flyjr ~ r gp as in 
IZhao et all l|2009l) than to r2oo as in lNavarro et al.l l|2004l ). 

But, at small r, the theoretical profiles become pro- 
gressively shallower and increasingly deviate from the NFW 
form (with constant inner logarithmic slope equal to — f). 
This behaviour is also in agreement with the results of nu- 
merical simulations showing that the Einasto profile gives 
a slightly better fit to the spherically averaged density pro- 
file for simulated haloes down to radii of about 0.00fi?vir. 
In Figure [5l we plot the Einasto profiles that better fit the 
theoretical profiles down to 0.01-Rvir. As can be seen, these 
Einasto profiles are still reasonably close to the theoretical 
profiles at radii four orders of magnitude less (r ~ 10~^i?vir). 

In Figure |6l the theoretical density profiles are also com- 
pared to the Eina sto profiles with m ass-dependent param- 
eters according to iGao et al] (|2008l ) and the same cosmol- 
ogy used by these authors (the mass-dependent Einasto pa- 
rameters are not known for other cosmologies). Although 
iGao et al] fitted the density profile of simulated haloes only 
down to 0.05 the total radiua^^l. the agreement between each 
couple of curves is similarly good as for the NFW profiles. 
In this case, the agreement of the theoretical profiles with 



iGao et al. I ll2008h used r200- To transform to the same radius 
i?vir as in Figure [5] their Einasto profiles have been extended to 
rgo where we have computed the halo mass. 
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Figure 5. Spherically averaged density profiles predicted by the 
present model (solid lines) for the same current masses M at iJvir 
as in Figures [2] and [3] (same colours, ordered as in Fig. [2]l, com- 
pared to th e dens ity profiles of the NFW form (dashed lines) with 
IZhao et al] ||2009|) concentrations fitting the average density pro- 
file of simulated haloes, both in the radial range covered by those 
authors (left) and down to a radius four orders of magnitude less 
(right) . For comparison we also plot the fits to the predicted den- 
sity profiles by an Einasto law (dotted lines). To avoid crowding, 
the curves for 10 Af» and 0.1 Af» have been respectively shifted 
upwards and downwards by a factor of 3. The lower panels show 
the residuals of the theoretical curves from the NFW profiles (left) 
and of the Einasto profiles from the theoretical curves (right) for 
the three halo masses (same order from top to bottom). 



the Einasto profiles at very small radii is not so good (al- 
though still much better than with the NFW profiles). But 
this is likely due to the rather limited radial range cov- 
ered in iGao et al] study. As the NFW profile gives, down 
to ~ O.OlJJvir, similarly good fits as the Einasto profile to 
the density profile of simulated haloes, the results shown in 
Figure [5] strongly suggest that the Einasto profile can do 
much better at small radii if the Einasto parameters are 
adjusted by covering a wider radial range. 

Given a protohalo with inner asymptotic density profile 

(3+q) ■ 



p oc Tp, so that Af(rp) oc Vp and £p{ 



(5 + a) 

OC Tp we 

have, from equation H24|) . that {p}{r) oc r". Since the pro- 
tohalo has a null logarithmic slope a (see Sec. [4]), it follows 
that the density profile for haloes must also have null inner 
logarithmic slope. In other words, there is strictly no cusp 
in the spherically averaged density profile for CDM haloes 
according to the present model. However, the finite central 
value is approached very slowly, in fact slightly more slowly 
than the Einasto profile (see Fig. [5|. 

The main differences th e density profiles p redict ed by 
the present model and by iDel Popolo et all (|2000t ) and 
lAscasibar et al.. ( 2004 ) models, both of which also use the 
SI framework, are that we do not assume spherical seeds 



Figure 6. Same as Figure [5] (including colours and order) but 
for the comparison with the de nsity profiles of the Einasto form 
(dashed lines) with [Cao et al.l (2008) mass-dependent parame- 
ters fitting the average density profile of simulated haloes in the 
ACDM cosmolog y adopted by those authors, both in the radial 
range covered by iGao et aU although extended out to iJvir ~ '"90 
(left) and down to the same small radii as in the right panel of 
Figure [5] (right) . The curves for 10 Mt and 0.1 M* have been 
shifted as in Figure [5] The lower panels show the residuals (same 
lines) for the three halo masses (same order from top to bottom). 



with convolved BBKS profiles and we make use of inside-out 
growth instead of an adiabatic invariant to determine {p){r). 
The model presented here is essentially equivalent to the 
Salvador-Sole ct al. ( 20071) model. The key difference is that 
Salvador-Sole et al.i (2007) used the typical cosmological ac- 
cretion rate onto haloes, whereas here we explicitly make 
use of the typical density profile of halo seeds. However, 
the typical halo mass accretion rate arises, as mentioned, 
from the typical protohalo density profile, so this difference 
is just a matter of presentation, motivated by the distinct 
theoretical framework of the two models: the peak and the 
excursion set formalisms. The only form al difference between 
the two models, apart from the fact that lSalvador-Sole et al.l 
also assumed spherical symmetry, is that the radius encom- 
passing a given mass adopted by these authors was inferred 
from equation ([24} although not using the spherical total 
energy of the re al protoobject but of its top-hat approxima- 
tion according to lBrvan fc NormanI (Il998l ) prescription. This 
should introduce additional numerical differences between 
the density profiles predicted by the two models. It is also 



important to mention that t he mo dels by IDel Popolo et al 



2000) and lAscasibar et all (|2004h and ISalvador-Sole et al 



(2007) include free parameters to be adjusted (the density 
contrast and filtering radius of the peak in the two former 
cases, and the value of AM/M setting the frontier between 
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Figure 7. Spherically averaged density profile for a halo with 
current total mass inside i^ yir equal to Mt predicted by the 
present model (g reen line), [ Salvador- Sole et al] | |2007| ) model 
(or ange line), the lAscasibar ^t al. ( 200 J" model (red line) and 
the iDel Popolo et aL I 1I2OOOI) model (blue line), respectively or- 
dered from bott om to top at the right end of each panel. The 
lAscasibar et al.l 's profile is only shown at radii larger than O.liJvir 
as recovered by those authors. 



minor and major mergers in the latter), while the present 
model includes no free parameter at all. 

In Figure [71 the spherically averaged halo density pro- 
file predicted by the present model is compared to those 
predic ted b y the [ Del Popolo et al.l ||2000| ). lAscasibar et al.l 
l|2004) and ISalvador-Sole et al.1 (|2007l ) models, all of them 
in the concordance model here considered. As expected, 
the theoretical profile predicte d by the present m odel is 
quite similar to that predicted bv lSalvador-Sole et al.l model, 
while these two profiles substantially deviate from those pre- 
dicted by the remaining models. On the other hand, the two 
form er profiles ar e in better agreement with the NFW profile 
with IZhao et al.1 ((20091 ) mass-concentration relation. More- 
over, the theoretical profile derived here is the one in best 
agreement with such a NFW profile, which is remarkable 
because it is also the only theoretical profile with no free 
parameter to be adjusted. 



6 IMPLICATIONS FOR MAJOR MERGERS 

The present model has been built for haloes formed by PA. 
Not only is this kind of halo formation crucial for the inside- 
out growth condition, but also for the possibility to apply 
the peak formalism in order to derive the peak trajectory 
leading to a typical halo with a given mass at a given time. 
The peak formalism itself is based on the existence of a one- 
to-one correspondence between haloes and peaks inspired 
in the SI model that ignores major mergers. This therefore 



raises two important questions. How do major mergers afi^ect 
the typical spherically averaged density profile for haloes 
derived under these conditions? And how do they affect the 
peak formalism? 

As shown in Section |3] given a seed with known spher- 
ically averaged density profile, we can find the spherically 
averaged density profile of the virialised halo evolving from 
it by PA, but the converse is also true. Given a halo grown by 
PA, we can calculate from its mass profile M(r) the spherical 
total energy of the protohalo, £p{M) (eq. [24]), and then de- 
termine its spherically averaged density profile, {pp){r), from 
equations (|29p and pOjl . Therefore, there is in PA a one-to- 
one correspondence between the initial and final spherically 
averaged density profil es. This is a well-k nown characteris- 
tic of spherical SI (Del Popolo et al.ll200d ). extended in the 
present paper to non-spherical SI. 

Interestingly, the reconstruction of the spherically av- 
eraged density profile for the seed of a halo having grown 
by PA can also be applied to a halo having suffered major 
mergers. This yields the spherically averaged density profile, 
(pp)(rp), of a putative peak that would evolve by PA into a 
halo with a spherically averaged density profile identical, by 
construction, to that of the original halo. Clearly, if the halo 
has grown by PA, such a putative peak exists and it is an 
ordinary peak. But if the halo has undergone major mergers, 
does it exist? Is it an ordinary peak? In which halo does it 
evolve? To answer these questions we will make use of the 
rigorous treatment of the peak formalism given in MSSa and 
MSSb. 

As mentioned, the peak Ansatz at the base of the 
peak formalism states that there is a one-to-one correspon- 
dence between haloes with M at t and peaks in the fil- 
tered density field at some small enough cosmic time ti, 
for some monotonous decreasing and increasing functions 
S{t) at Rf{M) of the respective arguments. According to 
this Ansatz, peaks associated with accreting haloes describe 
continuous trajectories in the 5pk-/?f diagram. Those peaks 
need not necessarily be anchored to points with fixed co- 
ordinates; they can move as the filtering scale varies (as 
real haloes do in the clustering process). But, thanks to the 
mandatory use of the Gaussian windowj^, the connection 
can be made in a simple consistent way between peaks trac- 
ing one given accreting halo at contiguous scale^ij. 

In a major merger, the continuous trajectories of (con- 
nected) peaks tracing the merging haloes are interrupted, 
while one new continuous peak trajectory appears tracing 
the halo resulting from the merger, leaving a finite gap in 
Rf due to the skip in halo mass between the merging and fi- 
nal objects. This is the only process where peak trajectories 
are interrupted. Haloes that do not merge but are accreted 
by more massive haloes are traced by peaks that do not 
disappear but become nested into the collapsing cloud of 
larger scale peaks tracing the accreting haloes. This leads to 
a complex nesting of peaks with identical 5pk but different 
R{. Once such a nesting is corrected, the number density of 



The decreasing behaviour of S with increasing R[ implied by 
the growth in time of halo masses is only guaranteed for that 
particular filtering window (see MSSa). 

A peak at scale R{ is connected with another peak at scale 
-Rf + dRf provided only these two peaks are at a distance smaller 
than R[ from each other (see MSSa). 
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peaks with (5pk at scales between R{ and Ri + dR{ and its 
filtering evolution recovers the mass function (and growth 
rates) of virialised haloes (MSSa and MSSb). 

Thanks to these results it can be shown (see MSSa) 
that the (5pk(t) and Rt{M) relations defining the one-to- 
one correspondence between (non-nested) peaks and (non- 
nested) haloes stated in the peak Ansatz are necessarily of 
the form (|37p - p8|l . These relations can be seen as the gener- 
alisation of those of the same form found in top-hat spher- 
ical collapse, with Sc{t) and q respectively equal to 1.686 
and 1. This does not mean, of course, that these parame- 
ters must take the same values in the peak formalism. On 
the contrary, the freedom in 5c{t) and q makes it possible 
to account for the change in the filtering window (Gaus- 
sian instead of top-hat) and possibly also in the departure 
from spherical collapse in the real clustering proces^3 For 
5c[t{z)] = 1.93 + (5.92- 0.4722 -h0.0546z2)/(l + 0.0005682^) 
and q — 2.75, the halo mass function predicted in the ACDM 
concordance cosmology (after correction for nesting accord- 
ing to MSSa) recovers the mass function derived from the 
excursion set formalism from 2 = up to any arbitrarily 
large 43- As the halo mass function at t predicted in the 
peak formalism is nothing but the filtering radius distri- 
bution for peaks with 5pk(t) (eq. [37]) transformed into the 
former by means of the relation (|38p . this result implies that 
there is indeed a one-to-one correspondence between peaks 
and haloes as stated by the peak Ansatz. 

The putative peak of a halo and its associated peak ac- 
cording to the peak Ansatz have the same 5pk(f) and R{{M). 
Their filtering evolution may be different: the trajectory in 
the 5pk--Rf diagram of the putative seed can always be traced 
down to any arbitrarily small filtering radius, whereas the 
trajectory of the associated peak can only be traced until 
reaching the filtering radius corresponding to the last major 
merger. However, the values Jpk and Rf of a peak do not al- 
low one to tell its 'past' filtering evolution because there is, in 
the filtering process, a kind of 'memory loss'. The Gaussian 
window, mandatory as mentioned for the peak formalism, 
yields a strong correlation between very close scales, which 
is welcome to carry out the connection between peaks at 
contiguous scales. But, at the same time, it yields a loss of 
correlation between scales different enough to encompass the 
gap produced in major mergers. Thus, owing to this partic- 
ular window, peaks at a given scale do not know whether or 
not they have appeared in some major merger (at a smaller 
scale). Therefore, the putative seed of a halo coincides with 
its associated peak according to the peak Ansatz and, as 
such, it is an ordinary peak. (It contributes to the filtering 
scale distribution of peaks with a given (5pk, regardless of the 
past history of the halo.) 

For this latter aspect to be true, we should adjust the mass 
function obtained form simulations or from the excursion set for- 
malism with values of Sc and q better adjusting the former mass 
function. In the present paper we fit, however, the mass function 
predicted in the excursion set formalism with Sc{t) = 1.686 and 
g = 1, so the departure from non-spherical collapse cannot be 
accounted for. 

The mass function predicted by the excursion set formalism is 
not fully accurate, particularly at large z, so one should rather fit 
the mass function drawn from numerical simulations. But this is 
not important here. 



Thus, according to the present model supported by the 
results shown in Section O the spherically averaged density 
profile of a halo having undergone major mergers must be 
indistinguishable from that arising from the evolution by 
PA of an ordinary peak (the putative peak) or, equivalently, 
there must be in the halo aggregation process a memory 
loss similar to that mentioned above affecting the filtering 
process of peaks. Such an implication has to do with the 
fundamental debate on whether or not virialisation is a real 
relaxation (e.g. iHen rikscn 20091). 

As noticed bv fPel Popolo et al.l l|2000l '). the one-to-one 
mapping between the initial and final density profiles in PA 
seems to indicate that there is no memory loss, at least in 
PA, during virialisation. This is at odds, however, with the 
fact that virialisation, even in PA, sets a time arrow. Even 
though the equations of motion of individual particles are 
time-reversible, owing to the highly non-linear dynamics of 
shell-crossing, any infinitesimal inaccuracies in the positions 
and velocities of particles in a simulation are rapidly ampli- 
fied, causing simulated orbits to deviate from the true ones. 
When the simulation is run forwards, this goes unnoticed 
because the system always reaches the same (statistically 
indistinguishable) final equilibrium configuration. However, 
when the simulation is run backwards it is readily detected: 
before shells can reach the maximum size (at turnaround), 
they begin to contract again towards the final equilibrium 
state. As the existence of a time arrow (or the impossibility 
of reversing by numerical means the evolution of a system) is 
an unambiguous signature of relaxation, we must conclude 
that virialisation is a real relaxation. 

In the case of PA the memory loss on the initial condi- 
tions is not complete: the order of shell apocentres is con- 
served. Consequently, despite the fact that information re- 
garding individual orbits is lost, we can recover the initial 
configuration in a statistical sense (in that we can recover 
the density profile). Thus, even in this case, the initial con- 
figuration cannot be exactly recovered because of the loss 
of information on the phase of particle orbits (individual or- 
bits are mixed up). And, in the general case including major 
mergers, the relaxation nature of virialisation is even more 
evident. The initial configuration of a given virialised ob- 
ject can never be recovered, not even in a statistical sense, 
because the density profile for virialised haloes does not har- 
bour information on their past history, so we can never be 
sure that the unconvolved density profile of their putative 
seeds, derived assuming PA, describes the density profile of 
their real seeds. 

The prediction made by the present model that the 
spherically averaged density profile of haloes does not har- 
bour information on their past aggregation history is thus 
consistent with the general behaviour of virialisation. More- 
over, it is consistent with several specific results regarding 
the density profile of simulated virialised haloes: i) virialised 
haloes with very different aggregation histo ries have similar 
spherically averaged density profiles (e.g. iNusser fc ShethI 
L1999:), ii) their density profile is ind ependent of the time 
they suffered their last major merger (|Wechsler et al.|[2002f ) 
and iii) their density profile does not allow one to tell how 
many, when and how intense the maj or mergers they have 
suffered are (|Romano-Dfaz et al.ll200g ). Certainly, there are 
also some claims in the literature that the density profile 
for haloes depends on their formation time. But this is due 
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to the particular definition of halo formation time adopted 
in those works, related to the form of the profile. Thus, the 
prediction of the model regarding the lack of information 
in the density profile of a halo on its past history is not 
contradicted by the results of numerical simulations. 

Thus, the present model strongly suggests that haloes 
having suffered major mergers have density profiles which 
are indistinguishable from those that have not. This predic- 
tion is important not only for the validity of the model (i.e. 
the validity of the inside-out growth assumption and the 
peak Ansatz), already supported by the good behaviour of 
the density profiles it predicts (Sec.O, but also because it 
allows one to understand why major mergers do not affect 
neither the typical halo density profile nor the peak formal- 
ism for halo statistics based on SI. Indeed, the fact that 
haloes formed in a major merger have, after complete virial- 
isation, would have a density profile indistinguishable from 
that of a halo formed by PA justifies that any halo with M 
at t can be associated with a peak with 5pk and Rf as if it 
had evolved from it by PA, as assumed in the peak Ansatz 
and that the typical halo density profile derived assuming 
PA is not affected by major mergers. As a corollary we have 
that the present model of spherically averaged density pro- 
file should hold for all haloes, regardless of their aggregation 
history. 



7 SUMMARY 

The present model relies on two assumptions: i) that the typ- 
ical density profile for virialised objects in hierarchical cos- 
mologies with coUisionless dissipationless dark matter can 
be derived as if all the objects grew by PA and ii) that the 
peak formalism correctly describes halo statistics. The key 
points of the derivation are that, contrarily to the kinematic 
profiles, the spherically averaged density profile does not de- 
pend on the (triaxial) shape of the object and that, during 
accretion, virialised objects grow from the inside out, keep- 
ing the inner structure essentially unchanged. 

Triaxial coUisionless dissipationless systems undergoing 
PA virialise by transferring energy from inner to outer shells 
through shell-crossing. Due to this energy loss, particles pro- 
gressively reaching turnaround describe orbits that contract 
orderly in the sense that the ellipsoidal surfaces effectively 
traced by their apocentres shrink and change their axial ra- 
tios without crossing each other until they stabilise. This 
causes the central steady object to develop inside-out. Such 
an evolution of accreting objects is in full agreement with 
the results of numerical simulations. 

One important consequence of the inside-out growth is 
that the radius encompassing a given mass in triaxial viri- 
alised objects formed by PA exactly coincides with the usual 
estimate of such a radius from the energy of the sphere 
with identical mass at turnaround. In this conditions, there 
is a one-to-one mapping between the density profile of the 
virialised object and its seed. This allows to infer the typi- 
cal spherically averaged density profile, {p)(r), of virialised 
objects from the energy distribution of their spherically 
averaged seeds, which can be calculated from the power- 
spectrum of density perturbations making use of the peak 
formalism. 

The consistency between the peak formalism and the 



one-to-one mapping between the initial and final density 
profiles following from PA implies that virialisation must 
be a real relaxation. This conclusion agrees with the general 
behaviour of virialisation (it sets a time arrow) as well as 
with specific results regarding the density profile of simu- 
lated haloes. As such, the spherically averaged density pro- 
file for haloes cannot harbour information on their past ag- 
gregation history which explains why major mergers do not 
alter the typical density profile of virialised objects derived 
under the PA assumption and do not invalidate the peak 
Ansatz. In other words, the sole condition that virialisation 
is a real relaxation is enough for the two assumptions of the 
model to be not only consistent but also fully justified. 

The model has been applied to CDM haloes. We have 
derived the typical unconvolved spherically averaged den- 
sity profile for peaks evolving by PA and from it the typical 
spherically averaged density profile for haloes with a given 
mass in a given epoch. Specifically, we have established the 
link between the typical halo density profile and the power- 
spectrum of density perturbations in any given hierarchical 
cosmology. The typical halo density profile so predicted in 
the ACDM concordance cosmology is in very good agree- 
ment with the NFW and the Einasto profiles fitting the 
spherically averaged density profile of simulated haloes down 
to one hundredth the total radius. However, such a theo- 
retical profile does not have a central cusp. In the ACDM 
cosmology, the model predicts that simulations reaching in- 
creasingly higher resolutions will find a typical halo density 
profile that tends to a core (null asymptotic slope) which is 
slowly approached, even slower than current fits using the 
Einasto profile. 
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